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Abstract. Let F(x,y) be a function of two variables, and suppose y = f(x) satisfies F(x,y) = for x in 
some range. L. Comtet and M. Fiolet ([3]) stated a combinatorial formula for the derivatives d n y/dx n in 
terms of the partial derivatives of F. Their formula, however, contains errors. Given the basic nature of this 
problem, it seems worthwhile to give the corrected formula, which is in fact slightly simpler than as stated 
in [3]- We give two derivations of the formula, by Lagrange inversion and by induction. We also correct 
Comtet and Fiolet's expression in [3] for the number of terms in the formula. 



1. Introduction 

Let F(x,y) be a smooth function of two variables, and suppose y = f(x) is a smooth function satisfying 
F(x,y) = for x in some range. It follows that dF — F x dx + F y ^ = where F Xl F y are the partial 
derivatives of F with respect to x and y. Hence 

(1) ^li^-El 

dx Fy ' 

as is well known. We can calculate the higher derivatives d n y /dx n ,n > 2 by taking the total derivative of 
equation (JlJ, giving 

d n y ,3 dy d d n ~ 1 y d F x d x F x 
^' dx~" = ^dx + fady^dx™- 1 ' ~ ^dx ~ Fy~dy* ~^ 

For any given n, this can be expanded to give an expression, analogous to equation ([1]), for d n y/dx n in terms 
of partial derivatives of F, but it is not clear how to give the resulting expressions as a uniform formula, valid 
for all n. L. Comtet and M. Fiolet stated such a formula in 1974 ([3J Theoreme 1], or see equations (|5Tj) -(j52" j) 
below), but their formula is not correct. Comtet and Fiolet indicated a proof that is undoubtedly correct in 
principle, using Lagrange inversion, but did not give their detailed calculations, so we do not know how the 
errors arose. 

The purpose of this note to give the correct result, which we do at equation ([Jj below. We also give two 
derivations, one using a form of Lagrange inversion (see equation (|12j) below), and another, elementary proof 
using induction. 

This note is organised as follows. In Section [2l we give some notation and in Section [3j we give our 
corrected version of Comtet and Fiolet's formula. A derivation using a Lagrange inversion method is in 
Section [4j and another by induction is given in Section In Section HI we discuss Comtet and Fiolet's 
original formula. In Section [Jj we observe that Comtet and Fiolet's expression [3j Theoreme 2] for the 
number of terms in their formula, is also incorrect, and we give the correct expression. 



2. Notation: Two-dimensional partitions 

To state our version of Comtet and Fiolet's formula, we need to introduce the concepts of one- and two- 
dimensional partitions, and set up some notation. A partition of the positive integer n, is a sequence of 
positive integers p± > P2 > ■■■ > p r > such that y^- pj — n. We write p h n to indicate that p — (p^ is a 
partition of n. The pi are called the parts of the partition, and the number r of parts is denoted \p\ . 

Analogously, if m and n are non-negative integers, not both zero, then a partition of (m, n) is a lexico- 
graphically ordered sequence 

(3) iplAl) > (P2, q2)~- > (Pr, Qr) > 0, 

where for each i, pi and qi are non-negative integers and not both zero, and JZiPi = m an< ^ = n - 

The lexicographic ordering is given by the rule (p,q) > (p',q') if p > p' or p — p' and q > q' . We write 
p h (to, n) to denote that p is a partition of (m, n) in the above sense. Confusion will not arise, since it is 

l 
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clear from what comes after the h, whether a one- or two- dimensional partition is indicated. Also as in the 
one-dimensional case, if p h (m, n) then we write \p\ to denote the number r of terms in the sum. 

For an ordinary partition p = (p,) h n, let e p ± denote the number of times the number k appears among 
the pi. Thus e p .k = \{i | Pi — k}\ . Naturally, for p h (m, n), we set e p ,fc,z = \{i \ pn = k and = l}\ , where 
Pi = (pn,Pi2)- That is, ep t fc j is the number of times (k, Z) appears in the partition. 

Next, we will encounter certain numbers a p where for an ordinary partition p h n, a p is defined by 

(4) 



p ILftilWi' 

and in the two dimensional case p h (m, n), by an obvious analogy, 
(5) 



where again the i term of p is being written as ijpiXiPa)- The coefficients a p in equation ([4]) appear in Faa 
di Bruno's well known formula, which we will use in our derivation (see equation 1 13[) . while those in equation 
([5]) appear in our main result at equation ([7]). ft seems noteworthy that the coefficients of these two formulae 
are analogous in this way. Finally, for a function g of one or two variables, we use the convenient notation, 

(6) 9 ^ = 11 dx^ ^ to = 11 dxPn dyPi2 - 

Again, note that it is clear from the context, whether p is an ordinary or a two dimensional partition. We 
continue to use g x and g y to denote the ordinary partial derivatives dg/dx and dg/dy respectively. 

3. The formula 

We now state our corrected version of Comtet's and Fiolet's formula. Let F(x, y) be a smooth function of 
two variables, and suppose y = f(x) is a smooth function satisfying F(x,y) — for x in some range. Then 
at a point where F y ^ 0, we have 

(7) 




The notation is as set out in Section 2; in particular, the coefficients a p are given by equation ([5]). See Section 
[5] for a discussion of where Comtet and Fiolet's formula differs from equation 0. The sum in equation ([7} 
is over two-dimensional partitions p h (n,r — 1) where \p\ — r, and (0, 1) cannot be a part of the partition. 
For example, taking n = 2, by inspection there are three partitions p h (2, \p\ — 1) with (0, 1) ^ p, namely 

(8) (2, 0), (1, f ) + (1,0) and (f , 0) + (1, 0) + (0, 2), 

and these have a p — 1 , 2 and 1 respectively, so we obtain 

r] 2 n F F F F 2 F 

(n\ " - xx 4_ 9 x x v _ r x r yy 

[ ' dx*- F y +A F% F* ' 

as may easily be verified using equation P)|. 

4. Derivation of equation ([7]) 

Comtet and Fiolet's indicated proof of their equation refers to an earlier paper [2], in which a further 
indication is given of a proof using Lagrange inversion. However, it seems more direct to use the following 
observation. In equation {J), provided F y ^ 0, F x /F y is the residue of F x /F at its pole in the y-plane at 
y = f(x). Therefore, 

(10) $ = _R eBv _ /(a) (^) = — L / d y^ 

v ; dx v_JWl F 2niJ c y F 

where C is a small loop around the pole at ?/ = f(x). The location of the pole varies with x, but for sufficiently 
small variation of a;, it remains inside C, so we may differentiate under the integral sign without moving C. 
Hence 

jn -i r an-l jp r)™- 1 F 

("» I?-5j/ *»^<t»- b -(^t<t»- 
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where the residue is at y = f(x) in the y-plane. We write this concisely as 

(12) ? = -Res,(^), 

where the subscript y refers to the residue being in the y plane. 

Equation (| 1 2j) is a relative of Lagrange inversion. The classic Lagrange inversion formula can be recovered 
by letting F[x,y) = x — g(y), so that / and g are inverse, and evaluating the residue. However, equation 
(fT2| is particularly suited to our present problem; all we have to do is evaluate the residue on the right hand 
side. To do this, we first introduce Faa di Bruno's formula. If z is a sufficiently differentiable function of y, 
and y is a sufficiently differentiable function of x, then Faa di Bruno's formula states that 

d n z ^ dMz 

( 13 ) ^=Z^^H^> 

phn " 

where the coefficients a p are given by equation |4]). The formula can be proved easily by induction. See p] 
Section 24.1.2(C)] or [5j for details and background. 

We can now begin our derivation of equation J?}- Applying Faa di Bruno's formula (| 13|) to logF(x, y) as 
a function of x, we have 

phn 

Hence equation (p~2|) becomes 

(is) ^ = EM~i) H (W-i)!iW^)- 

phn 

Now to simplify the notation, we temporarily work at the origin x = y = 0, assuming therefore that 
F(0,0) = and F y (0) ^ 0. Then the pole of F x>p /F^ at y = has order \p\ , so 

(16) Res^(f^) = ilpl l iy jll^ Am- W )\ v =o. 

By Leibniz' rule, 

flbl-i '^Iz 1 Clnl-lV ^H-i-sp ,9 s 

< 17 > 5H=rW-W»)- M ) ^ g .k^-!-.), ^---' ^"^ '■ 

and by Faa di Bruno's formula (fl3|) applied to (F/y)~^ as a function of F/y, we have 

(is) JL(( F / v) -\P\)-x^ a (.pki Qpl + M-i) 1 (g/»)M 

U8j dy° (( /V) >~^ s q[ j (H-l)! 

To take the limit as y tends to 0, in equation (118|) . we replace d r (F/y)/dy r with ;qrj-<9 r+1 F/dz/ r+1 for r > 0. 
Hence F/y becomes F y and 



(19) (F/y) v , q = 11(^7^) becomes Y\ t>l { qi + l) [Fy)y ' q 



i>X 

Therefore at y = 0, equation (fT%j) becomes 



— ((F/,)-H)- V« M) M (IPl + M-1)' (^)»,« 



" ^ "■' ' (I/) - 1 )! I !•((/; - I i T7\l'\ ' ■< 

qhs 

Substituting this into equation (|17j) and simplifying gives 

(m Hrnr ^'" 1 r F «* ^ = V " 9 (-i) M (|p| + |g|-i)! (jkki^T^ 

1 j B-VfylPl-i^/yjW" ^ s!nfe + l)(b|-l-s)!4 p|+l91 dyH- 1 -* ' 

Then by equations (fT5|) and (fT6|) . 

^ = V (-l)H+l^lq p a g (|p| + |gj - 1)! (F y ) y , q d^^F^ 

[ ' dxn pVn s\Y\{ qi + l){\p\-l- s)\ jpM+lsl ' 

ghs< |p| — 1 
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Although we were assuming x = y = in the derivation of equation (122p , the same equation clearly holds at 
any point. We use Leibniz' rule to write out 

d^-'F^ _ dH- 1 - n d Pl F _ v (| p | _ i _ 3 )j Igl ffli+ftf 

1 ' flwlpl-i-* ; ^ /3i!.../3i„| ! H flaP^A ' 



Then equation (|22)l becomes 

(24) - V (-l) |p|+l9l a P a g (b| + |gj - 1)! ( F ak<? 

ghs< |p| — 1 
l 8i,..., j 9 b |:Eft=|p|-l-» 

In equation (|24[) . each term corresponds to a partition p h n, a partition g h s for some s < \p\ — 1, 
and \p\ non-negative integers /3i, 0\ p \ with ^ /3j = \p\ — 1 — s. For a given p, q and (/3j), there is a unique 
permutation 7r of 1, .., \p\ such that the resulting sequence (p^-m, /3 n n)), (p-n-n p \), Arflpl)) is lexicographically 
ordered, and p, /3 and 9 thus can be combined in a unique way into a two-dimensional partition P(p, /3, q) 
given by 

lexicographically ordered 



(25) P(p,(3,q) : (p^i), Ar(i)), (?MH), &t(H))j (0,<7i + 1), (0,q\ q \ + 1)- 

We see that P = P(p, /?, g) is lexicographically ordered throughout; indeed p^(\ p \) > 0, so (p v (\ p \), /^(Ipl)) > 
(0, Qi + 1). Also P satisfies |P| = |p| + |g| , P h (n, |P| — 1) and (0, 1) ^ P. Furthermore, clearly any partition 
satisfying these two relations, is of the form ([23]) . Moreover, 

r „ fi s ^Pi jy ft--F 1 .P| P |/ b |(^) j y,g _ P X ^ P 

f IpI+i?i - ~~^T' 

In particular, this term is the same for each triple (p, (f3i),q) in equation (|24[) which gives rise to P, so we 
can change the summation to be over partitions P, which will lead us to the final form required. To do this, 
we need to calculate how many triples (p, q) give rise to the same P. But P occurs as P = P(p, (/%), q) 
for a unique p and q, and for permutations of the (/%) which preserve the set 

(27) {Oi,A),-,(Pb|,^b|)}, 

so that after lexicographic ordering, the same partition P is obtained. The number of such permutations is 

(28) n^/ n 

fc fc>0,2>0 

(The product in the denominator is over k > 0, as the multiplicities ep$ t i are not relevant, being those 

associated with the part of P coming from q.) All the associated terms in equation (|2~i|) are the same, so 
equation (|24|) becomes 



f291 d^y v (-l)^a P a q (\p\ + |gj - 1)! U k e p , fc ! P^p 

dx " ph(^i-i) /3i!-/3|pi^nife+i) n fe >o, ; >oe W ! pi p| ■ 

The coefficient in equation ([29]) is 

(-l)W+Ma pQg (|p| + |g|-l)!n fc e M ! _ (-1)H+Mn!«!(bl + | 9 | - 1)! n fc e M ! 



>0,i>0 e P,k,l- 

(30) = (-l)l p ln!(|P|-l)! 

IL Pii { - Hi P2! I!, e Pl0 j\ Yl k >o,i> ep.fc,;! 
= (-1) |P| «p, 

where we have used that Y[ Pi! = YiPil an d II = II A- life + !)■> which are evident from equation 
and the definitions of a p from equations ([?]) and ([5]). Thus equation (12T))) amounts to equation ([7]), and the 
derivation is complete. 
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5. Proof of equation (J71) by induction 

It is possible to prove equation ([7]) by induction on n. The inductive proof is completely elementary and 
it therefore seems worthwhile to give this proof here. To simplify the presentation, wc will find it useful to 
introduce the following operations on partitions. For a partition p b (n, m), and for i,j with i > 0, j > 0, 
(hj) 7^ (1 5 0) and e Pj ij > 0, define q = Aj 3 -(p) to be the partition with 

{e p ,fc,i-l (k,l) = (i,j) 
e P ,k,i + 1 (k,l) = (i - 1, j) . 
e p ,k,l otherwise 

In words, viewing the table of muliplicities for p, we move one unit from the position to the (i — 
position. Notice the condition ^ (1,0) ensures that the operation results in a partition, i.e. we have 
e Ai ,(p),0,o — 0. Similarly, for i, j with i > 0, j > 0, (i, j) ^ (0, 1) and e p ,i,j > 0, define q = Bij(j>) to be the 
partition with 

{e p ,k,i ~ 1 (k,l) = (i,j) 
e p ,k,i + l (M) = (U-l) • 
e p ,k,l otherwise 

The operations Aij and Bij are convenient for expressing the derivatives of a term F Xt y tP , via the following 
two formulae. 

^ 11 [ dx k dy l) ' ' ih-i),,- ; ' ' dx^dyi 



(33) = £ 



i,j:e 9 ,i_i,3>0 

= ^ ] ( e P ,i-l,j + l)-Fx,j/,p, 

P,»,j':Ai, 3 (p)=9 

where the sum is over all (p,i,j) such that _Ajj(p) is defined and equal to q. Similarly, 

d 

(34) -Q-(F x ,y.q) = ^ ] ( e p,i,j-l + l)^x,l/,p- 

Suppose p b (n, to). Then it is easy to verify that 

(35) IAi(p)| = |Bij(p)| = bl 
and 

(36) Ajj(p) b (n — 1,to) and Bij(p) b (n,TO — 1), 
and further, using the definition at equation ([5]), 

(37) " /'•'■.' .. i 



= ««,,., ;.; . i)"" and = m( e PM :; + D Qp - 

We also need the following, where the backslash notation indicates removing the specified parts from a 
partition p b (n, m) containing at least one of each part: 

/oo\ e p,i,o e p-i,i j 2e Pj i ) oe Pi o.2 

(38) a P \(i,o) = — ^- a p ,a p \ (M) = ~^a p and Op\(i,o),(o,2) = nm{jn _ 

Combining the second equation of (|37[) and the first equation of (|38|) gives, where again p b (n, to), 
Bi,3(p\(i,o)) TO(e p \( 1)0 ),i,i_i + 1) p ^ (1,0) TOn(e p \ (lj0)iij _ 1 + l) p ' 
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wherever Bi^(p\(l,0)) is defined. (In the right hand we have replaced e p \(i t o).i.j with e p .i.j, since they are 
equal for j > 0, and Bij is only defined for j > 0.) 

We are now ready to begin the proof of equation by induction. For n = 1, equation {7|| is just equation 
(1), so by induction assume equation {2} holds for n — 1. Then by equation ([2]), we have 

d n y 



(40) 



qh(n-l,| 9 |-l) 
(0,1)^9 



For any partition q in the sum, we have four terms 



(41) 
where 



d_ _ F±d_ 

dx F y dy 



KFy l9l F x>y<g ) =X 1 +X 2 + X 3 + X 4 



(42) 



X 1 

x 4 



dx { y ' x ' v ' 9 



f-kl-ij? 



F-\i\ — (F ) 

y dx [ x ' v ' q) 



= F~w 
y 



F y dy 



\q\F-M- 2 F x<y . 



\&p,i — l,j "b f)-^a:,y,pi 



gU(l,0)U(0,2)i 



y 



d 



(e p \(l,0),i,3-l + l)-^x,i/,p) 

p,t,j,B 1 ,j(p\(X,0))=g 



and we have used equations (|33|) and (|34|) respectively, to obtain the expressions for X 2 and X4. The sums 
in X 2 and X4 are thus over all triples (p,i,j) where p is a partition and i,j are integers, satisfying that 
A i,j{p) = Q respectively B itj (p\(l,0)) = q. 

In equation (|42|) , we see four types of relationship between partitions p and q, corresponding to the four 
terms X\-±. The following lemma collates information about partitions that are related in one of these ways. 

Lemma. Let p and q be two-dimensional partitions and let p b (n, m). Suppose that p and q are related by 
one of the following four equations. 

1. g = p\(l,l). 

2. q = Ai_j{p) for some 

3. g = p\(l,0),(0,2). 

4. q = Bij(p\(l,0)) for some 
Then the following hold: 

(1) We have 



\P\ = 




in case 1 
in case 2 
in case 3 
in case 4 



and q b 



(n — 1, m — 1) 
(n — 1, to) 
(71 — 1 , to — 2) 

(71 — 1, TO — 1) 



in case 1 
m case 2 
m case 3 
m case 4 



(2) p b (n, |p| — 1) i/ and only if q b (n — 1, g — 1). 

(3) // (0, 1) ^ g i/ien (0, 1) ^ p, and if (0, 1) ^ p then (0, 1) ^ g unless either p and q are related by case 2 
with = (1, 1), or p and g are related by case 4 = (0,2). In either of these exceptions, 
we always have (0, 1) G g. 

(4) 

in case 1 
in case 2 
in case 3 
m case 4 



nm(m-l) "P 

jep,i,je p ,i, 



.J-i 



+ 1) 



Proof. (1) The assertions about \p\ are clear in cases 1 and 3, while cases 2 and 4 are immediate from equations 
(|35| . Similarly, in cases 1 and 3 it is clear that g b (n — 1, m — 1) and g b (n — 1, to — 2) respectively, while 
equations 136[) show that ph (n — 1, to) and p b (n — 1, m — 1) in cases 2 and 4 respectively. This gives (1). 
(2) follows directly from (1). 
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(3) First suppose (0, 1) ^ q. Then it is clear in cases 1 and 3 that (0, 1) ^ p. For case 2, we need to show 
that (0, 1) € p implies (0, 1) € Aij(p) whenever the latter is defined. This follows from the definition of A+j 
at equation (f3"T|) . Similarly, for case 4, we need that (0, 1) 6 p implies (0, 1) £ Bij(p) whenever Bij(p) is 
defined, which again follows from the definition |32|) . 

Conversely, suppose (0, 1) ^ p. Then again it is obvious in cases 1 and 3 that (0, 1) ^ q. In cases 2 and 4, 
we have from the definitions (|31[) and (f3"2"| . that (0, 1) ^ p implies (0, 1) ^ Ajj-(p), unless (i, j) = (I, I), when 
(0, 1) € ^4i,j(p), and similarly (0, 1) ^ p implies (0, 1) ^ B it j(p), unless (i, j) = (0, 2), when (0, 1) £ B itj (p). 

(4) merely collates the second equation in (|38[) . the first equation in {37]), the third equation in (f38|) and 
the last equation in (|39|) . for cases 1 — 4 respectively. □ 



We now complete the proof by induction. First, by comparing the powers of F y appearing in each term 
X1-4 in equation (|41|) . with the assertions about \p\ in part (1) of the lemma, we see that each term is of 
the form F y F Xjy>p , where p and q are related by one of the cases 1 — 4 of the lemma. By parts (2) and 
(3) of the lemma, since q h (n — 1, \q\ — 1) and (0, 1) ^ q, we have p h (n, |p| — 1) and (0, 1) ^ p. We have 
therefore established that equation (|4^)) is of the form 

ph(n,|p|-l) r a 

for some coefficients i p . It remains to show that t p — (— l)' p 'a p . To do this, we fix a particular p h (n, m) 
where m = |j>| — 1 and (0, 1) p. By equations (|40| and (|4Tj) . we have 

(44) i p = A + B + C + £>, 

where A, i?, C, D are respectively the contributions to t p from all gh(n-l, \q\ — 1) with (0, 1) ^ g, via the 
terms Xi_4 in equation (|42]) . These are given by: 

A = -(-l)' 9i a g \q\ if (1, 1) e p and A = otherwise, 

B= (e^-^ + lX-l) 19 ^, 

q,i,j:A i , j (p)=q:(i :] )^(Q,2) 
45 ) I 

C* = (~l) lql a q \q\ if (1, 0), (0, 2) E p and C* = otherwise, 



D = - Y, (^(1,0)^-1 + l)(-l) kl OflM • 

g,i,j:B 4 ,j(p\(l,0))=g:(t,i)^(l,l) 

Note that in equations (|43|) . we have the conditions (1, 1) ^ p and (0, 2) ^ p respectively in the sums for B 
and D. By parts (2) and (3) of the lemma, these are the conditions needed to ensure that (0, 1) ^ q and 
q h (n - 1, |g| - 1). 

Using parts (2) and (4) of the lemma, we compute 

(46) A = -(-1)1*1 M ^ia p = (-l)H(bl - l)te a = (_i)H£Eii a 

nm n 

noting that this equation automatically gives the correct zero contribution when (1, 1) ^ p. Next, 

(47) B= £ + *y . op- £ (-D H ^f«P- 

Note that the range of summation in equation (|47|) is the set (i,j) having Ai t j(p) defined and ^ (1, 1). 
By equation (f3"Tj) . Ajj-fjj) is defined where i > 0,j > e p ,i,j > and (i,j) ^ (1,0). However if i = or 
e P ,i,j = then the corresponding summand is automatically zero, so the conditions i > and e p ,i,j > are 
removed in the right hand expression. 
Similarly to equation (|46p . we have 

(48) C = (-l)W 2e ;'° e °'> p |g| = (-ijW^o^W 

nmym — 1) nm 

and similarly to equation (|47|) . 



(49) ° = E -(e P \(i,o)„--i + l) nm( } ^| g |= E ^ 

(i,jV(l,0),(0,2) 
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Substituting in equation (|44[) gives 

^ _ (_i)Hq, ) e P' 1 > 1 _|_ ^ e p,j,j _|_ ^ e p,i,o e p,o,2 _|_ jep,i,pep,i,j 

p p ' n ' n nm ^— ' nm 

(i,i)^(l,0),(l,l) (*,j)^(0,2) 

= (-i)MoJ E ^+E^? M 



(-i) h «p{ 



(».i)^(l,0) i,j 



n n > 

= (-l) |p| a P , 

where we have used ^ i ■ ze Pi i j = n and . je p ,i,j — m. This completes the proof. 

6. Discussion of Comtet and Fiolet's formula 

In this section, we compare Comtet and Fiolet's formula in [3] with our equation ([7]). Comtet and Fiolct 
write E = N x N\ {(0, 0), (0, 1)} where N denotes the natural numbers including zero. Also, in a notation 
slightly different from ours above, they let F^j denote d i+:> ' F ' jdx % dyK Then their formula ([3] Theoreme 1]) 
is 

rl n ii 2 "~ 1 _1 

< 51 > 3?=E<^>-'..- 

m—l y 



where 



T _ \- n\q\(q)s{q + S) C1 -pr 

{ ' n ' m ~ 2^ FT TkV^+h 11 



h+2h+3h+-=n Alfe > lV ; (j,i)e-B J 

ci+2c2+3c 3 + ...=m — 1 

the sum being over tables of non-negative integers (dij) : 6 E, and where the following notation is 

used: 

/, = dij; Cj = di t j (row and column sums), 

i » 

(g)s = 5(5 + i)...(g + 5-l). 

To see the basic correspondence between equations (l52j) and (J7J, we regard each table of non-negative 
integers (dij) : € E in equation (|52p , as the multiplicities of a unique two-dimensional partition p, so 
that dij — e.p,i,j- We have (0, 1) ^ p, by definition of E. The conditions under the sum in equation (|52[) arc 
then equivalent to requiring that p h (n, m — 1). To assist comparison, we rewrite our equation ([7]), to look 
similar to equations (|51[) and (|52[) , putting m = \p\ to give 

(54) ^ = E(^r E 

m=l y ph(n,m — 1) 

|p|— 771 

The upper limit of summation m = 2n — 1 in equation l|54p is automatic, since if p \- (n, |p| — 1) and (0, 1) ^ p 
as required in our formula, then \p\ <2n — 1. To see this, write p in the form 

(55) (pi,Qi), ■-, (Pr,q r ), (0, q r +i), -, (0, q r +s), 

with p r > 0. Then \p\ = r + s, and since c^+i, q r + s > 2, we have r + s — 1 = \p\ — 1 > 2s, so r — 1 > s, 
while also n > r, so 2n — 1 > 2r — 1 > r + s = \p\ . 

Now for a partition p, to assume Comtet and Fiolet's notation, write dij instead of e P ,ij- Then 

(56) n^ i! = nw* and Hp^ [ = n^s 
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where Ck and Ik are defined in equation (|53[) . Thus since ph (n, m — 1), we may write equation (|54|) as 



2 v^ 1 / _ l\m \ - n!(m — 1)! -p-r 



1 ' dx n ^ K FJ ^ UAkl) Ck+lk 11 dal 

m=l " /i+2/ 2 +3i 3 + ...=n 1 lfc v ' ijeB 



Ci+2c2+3c 3 + ...=m — 1 

There are just two differences between our equation (|57[) and Comtet and Fiolet's formula comprised in (J5T 
and J5 



(1) We have the additional condition Iq + 1\ + .. = . djj = m under the inner sum. 

(2) We have a term (m — 1)! in the coefficient, where Comtet and Fiolet have the term q\(q)s{ ( l + S) Cl . 
We deal with these points in turn. (1) Our additional condition Iq + l\ + .. = ■ <iy = m under the 
inner sum, ensures that \p\ = to, or in other words that p h (n, \p\ — 1). There is no equivalent condition in 
equation (|52[) . but in the text immediately following their statement, Comtet and Fiolet observe that the 
polynomial I n ,m defined by equation (|52p is homogeneous of degree to, in the sense that ^ dij = to. This 
is not an automatic consequence of the other conditions: Indeed, equation (|52[) as it stands includes terms 
for any partition p h (n,m — 1) with to < 2n — 1, and these certainly do not all satisfy \p\ = m. (Perhaps 
the extra condition was simply lost during the typesetting of (|52|) ; there is a row of dots under the condition 
ci + 2c2 + 3c3 + ... = m — 1, in equation (|52[) as it appears in [3], which look like separators for a further line 
below that was not printed.) 

Turning to difference (2), Comtet and Fiolet have a term q\{q)s(<l + S) Cll where in our equation (|57[) we 
just have (to — f )!. However, these terms are nearly equal (but our version is correct, as we check below with 
an example). From the definitions at equation (|53[) we have 

(58) q + S + ci = 1 + ^ jc i+ i + Cj + ci = 1 + y^jCj = to, 

j>1 j>2 j>1 

and hence in equation (j52"f . 

(59) g!(«) s (g + S) C1 = q\(q + l)...(q + S)(q + S + l)...(q + S + c x - 1) = q(m - 1)!. 
Thus, equation (fB"2")) becomes 

(60) / - V n\q(m-\)\ yr <^ 

h+2l 2 +3l 3 + ...=n iife>lV ^ J 
ci+2c 2 +3c3 + ...=m— 1 
io+il+'2---=ni 

Now, the terms of equations (|60|) and (|5T[) agree except that there is an additional factor of q in (|60f . To 
check that our version (|57[) is indeed the correct version, consider the case n = 5, and the coefficient of the 
term belonging to the partition p h (5, 4) given by 

(61) (5, 4) = (1, f ) + (f , f ) + (1, f ) + (f , 0) + (1, 0) + (0, 2). 

Note that this partition satisfies p h (5, \p\ — 1) and (0, 1) ^ p. The coefficient of the corresponding term 
F^F^ y Fyy in the expression for d 5 y/dx 5 is 600. This is stated correctly by Comtet and Fiolet in the numerical 
table in [3] (it is the coefficient of /f 1/0,2 m their notation), and we have confirmed its value by performing 
a recursive calculation up to n — 6 using equation However for this partition, we have 

(62) di,o = 2,di,i = 3,do,a = l, 

and all other dij are zero. Hence the terms going into Comtet and Fiolet's formula are as follows: 

n = 5, 

'0 = 1> '1 = 5, 1% = ... = 0, 

(63) co = 2,ci=3,C2 = l,c 3 = ... = 0, 

S = C2 + C3 + ... = 1, 

q = 1 + c 2 + 2c 3 + ... = 2. 

aso the coefficient according to Comtet and Fiolet's formula (I52| is 



5!2!(2) 1 (3) 3 _ f20x 2x2x3x4x5 _ ionn 
(b ' (2!)if!2!3! ~ 4x6 ~ ' 
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which is out by a factor of q = 2, as expected, while our formula © gives the coefficient correctly as 

5!4! 



(65) 



600. 



p 2!2!3! 

In summary, Comtet and Fiolet's formula, once corrected, should look like the equation below, which differs 
only notationally from equation l|7|. 

2n-l 



(66) 



dx n ^ 1 P 



1. 



E 



n!(?-l)!(g) s (? + S) cl 



m=l y h+2i2+3i 3 + ...=n 

Ci+2c2+3c3 + ...=m— 1 
co+ci+C2+...=m 



n fc >i(fc! 



I^Cfc+ifc 



F 

II 



(ij)£B 

with the proviso that the expression (q — l)\(q)s(q + S) Cl is just a complicated way of writing (m — 1)! 

7. The number of terms in equation (JTj) 

Comtet and Fiolet gave a formula ([3j Theoreme 2]) for the number of distinct terms appearing in their 
formula (|51fl52|) . as a function of n. In terms of our equation [[7]). this is the number of partitions p h (n, \p\ — 1) 
with (0, 1) ^ p. Let this number be a{n). Comtet and Fiolet gave the following table of values of a{n) for 
n < 23 (we have added a(24) from our own calculations). 



n 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


a(n) 


l 


3 


9 


24 


61 


145 


333 


732 


1565 


3247 


6583 


13047 


n 


13 


14 


15 


16 


17 


18 


19 


20 


21 


22 


23 


24 


a(n) 


25379 


48477 


91159 


168883 


308736 


557335 


994638 


1755909 


3068960 


5313318 


9118049 


15516710 



The sequence 1, 3, 9, 24, 61, ... of values a(n) can be found in the On-Line Encyclopedia of Integer Sequences 
[5J. The numbers a(n) were correctly stated by Comtet and Fiolet, but their formula Q3, Theoreme 2]) 
is not correct. This formula states that a(n) is the coefficient of ^"w" -1 in ri(i j)e.E(l ~ fu 1 )^ 1 , where 
E = N x N\ {(0, 0), (0, 1)} as defined above equation fBTj) . This is seen to be incorrect, for example, 
by calculating the coefficient of t 2 u. Expanding j)e£;(l — t 1 ^)' 1 , to all the terms that can possibly 
contribute, we get 

(67) (1 + t + t 2 ){\ + tu)(l + t 2 ){\ + t 2 u), 

from which we see that the coefficient of t 2 u in JT^ j)eE^ — t 1 ^)^ 1 is 2, while a(2) = 3 from the table above 
(or from equation @). Comtet and Fiolet's incorrect formula is repeated on page 175 of [3]. The correct 
result is as follows. 

Theorem (Compare with [3], Theoreme 2 or [J], pl75). The number a(n) of partitions p h (n, \p\ — 1) and 

(0, 1) ^ p, or equivalently, the number a(n) of distinct terms in equation Q) ; is given by: 



a(n) = Coefficient of fw"" 1 in TT - 

11 1 — t l u 



Proof. Write 

( 69 ) n*,«)= n !_ t ^ +J -i = n (H^ i+j - x ) d ) 

(i,j)eE (i,j)eE d>o 

For a given partition p h (n, |p| — 1) with (0, 1) ^ p, let <ii j:) - = e Pj jj ; £ S, be the table of multiplicities 
of p. Corresponding to this p, we pick the term with d — dij in the factor of F, and get a term of F 
equal to 

(70) ^. i , j idi,j v E, iij {i+j-l)d i , j _ 

This term contributes 1 to the coefficient of t n u n ~ x exactly when 

idij = n and ^~~^(i +j— = n — 1 
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(71) ^ id tj = 11 and ^ jdjj = ^ di 'i ~ 

that is, exactly when p h (n, \p\ — 1) as required. □ 



The numbers a(n) can be calculated easily from the generating function (|68[) . and we indicate the simple 
scheme we used to recalculate the figures shown in the table (which, as we mentioned above, agree with the 
values given by Comtet and Fiolet for n < 23). Write 

(72) F(u,t)= [] j— f^I = J2Pn(t)u n , 

(ij-)GE n>0 

so that a(n) is the coefficient of t n in p n _i(t). Also let 

(73) G(u,t) = logF(u,t)= -log(l - tV+^ 1 ) = Yl {tu%+ °~ 1)r = q m (t)u m . 

(i,j)eE (iJ)EE,r>0 ' rn>0 

For m > 0, we see that 

m/d+i .id 

(74) q m (t)= £ fm /{i+j -i )= £ J- '-j, 

-i,j>0 d divides rn i— 

— 1 divides m 

from which the q m (t) are easy to compute explicitly. Now F = e G , so differentiating we have -FG U — _F U , or 



(75) EM*)« r )(X>.(*)«* _1 ) = J^Pnit) 



nu n -\ 



r>0 s>0 n>0 

n-1 



and comparing the coefficient of u n gives 

1 - 

(76) R,(*) = -E 



SQsPn — s • 

n * — ' 

s=l 

Starting frompo(i) = t) = 1/(1 — t) = ^2 n > t n , equations GH) and (|76p enable successive calculation of 
the polynomials p n (t) for n > 0, and then a(n) is extracted as the coefficient of t n in p n ^i(t). If we desire the 
values of a(n) for 1 < n < E, say, then we only have to keep terms up to t E in the polynomials p and q. To 
allow the reader to verify that this calculation does produce the numbers a(n) in the table above (and hence 
also the numbers stated by Comtet and Fiolet), below is some self-explanatory VBA code implementing the 
calculation of a(n) for n < 50, say. 

Sub Calc_AofN_upto_EO 
E = 50 

ReDlm p(0 To E - 1 , To E) : ReDim q(0 To E - 1 , To E) 

For m = 1 To E - 1 : For d = 1 To m 

If m = d * IntCm / d) Then 

For i = To • / d t 1 

If d * i <= E Then q(m, i * d) = q(m, 1 * d) + 1 / d 
Next : End If : Next : Next 

For j = To E 
p(0, j) = 1 
Next 

For n = 1 To E - 1: For s = To n: For j = To E: For i = To j 

pCn, j) = p(n, j)+l/n*s* qCs, j - i) * pCn - s, i) 
Next : Next : Next : Next 

For n = 1 To E 

Debug. Print pCn - 1, n) 

Next 
End Sub 
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